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Abstract 

We present complete solution of Altarelli-Parisi (AP) evolution equation in next-to- 
leading order (NLO) and obtain t-evolution of non-singlet structure function at low-x. 
Results are compared with HERA low-x and low-Q^ data and also with those of com- 
plete solution in leading oder (LO) of AP evolution equation. 
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1. Introduction: 

The Altarelli-Parisi (AP) evolution equations [1-4] are fundamental tools to study the 
t(=ln((5^/A^)) and x-evolutions of structure functions, where x and are Bjorken scaling 
variable and four momenta transfer in a deep inelastic scattering (DIS) process [5] respec- 
tively and A is the QCD cut-off parameter. Though numerical solutions are available in the 
literature [6] , the explorations of the possibility of obtaining analytical solutions of AP evolu- 
tion equations are always interesting. In this connection, complete solutions of AP evolution 
equations at low-x in leading order(LO) have been obtained [7]. Its natural improvement 
will be the next-to-leading order (NLO) calculation. 

In this paper, we present complete solution of AP evolution equation in NLO at low-x and 
obtain t-evolution of non-singlet structure function. Results are compared with the HERA 
low-x low-Q^ data, and also with those of complete solution in LO. Here Section 1, Section 
2 and Section 3 give the introduction, the necessary theory and the results and discussion 
respectively. 



2. Theory: 

The AP evolution equation for non-singlet structure function in NLO is [8] 
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where, 



f{w) = fPpiw) + 2Pg{w) + InjP^.iw) + '^Pa{w). 



The exphcit forms of higher order kernels are [9] 
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Running couphng constant in higher order has the form [10,11] 
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being the number of flavours. 

Using Taylor expansion method [12] and neglecting higher order terms as discussed in 
our earlier works [7,13,14], F^^(x/w,t) can be approximated for low-x as 
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Putting equation (2) in equation (1) and performing u-integrations we get. 
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where, 



A(x) = -{-2xlnx + x{l - x^)}, 
o 

1 — w 

B{x) — ^ J f{'w)dw, 

C{x) = ^{3 + 4/n(l - x) + {x - l){x + 3)} 



and 



D(x) — — f(w)dw + x / f(w)dw. 
Jo Jo 

For a possible solution, we assume that 



2% 



k 



2% 



where, k is a numerical parameter to be obtained from the particular Q^- range under study. 
By a suitable choice of k we can reduce the error to a minimum. Now equation (3) can be 

recast as 
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The general solution of equation (4) is 

F{u, y) = 

where, F is an arbitrary function and 

U{x, t, Fi^^) = Ci and V{x, t, F^^^^) - C2, 
where Ci and C2 are constants, form a solution of equations 
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Solving equation (5) we obtain. 
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If U and V are two independent solutions of equation (4) and if a and/3 are arbitrary 
constants, then 

V = aU + p 

is called a complete solution of equation (4). Then the complete solution [12] 
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is a two-parameter family of planes. The one parameter family determined by taking (3 = 
has equation 



Differentiating equation (6) with respect to a, we get 
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Putting the value of a again in equation (6), we obtain 
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which gives the t-evolution of non-singlet structure function F^^{x, t) in NLO. In an 
earher communication [7], we suggested that for low-x in LO 
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We observe that if b tends to zero, then equation (8) tends to equation(9), i.e., solution of 
NLO equation goes to that of LO equation. Physically b tends to zero means number of 
flavours is high. 
Again defining. 
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which gives the x-evolution of non-singlet structure function F^^{x,t) in NLO. 

Proton and neutron structure functions measured in deep inelastic electro-production can 
be written in terms of singlet and non-singlet quark distribution functions as 

FUx,t)^^Fi{x,t) + ^Fi^'{x,t) 



and 



These equations give 
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from which we can calculate experimental values of F2 in t and x ranges given in and 
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3. Results and Discussion: 

In the present paper, we compare our results of t-evolution of non-singlet structure func- 
tions from equation (8) with the HERA low-x and low-Q^ data [15]. Here proton and 
neutron structure functions are measured in the range 2 < < 50 GeV'^. Moreover here 
Pt < 200 MeV, where Pt is the transverse momentum of the final state baryon. We con- 
sider number of flavours rif—A:. 

In figures l(a-c) we present our results of t-evolution of non-singlet structure functions 
F2^'^ (sohd hues) for the representative values of x given to test the evolution equation (8) in 
next-to-leading order. Agreement is found to be good. In the same figures we also plot the 
results of t-evolutions of non-singlet structure functions ^2^*^ (dashed hues) for the complete 
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solutions from equation (9) in leading order. Data points at lowest-Q^ values in the figures 
are taken as inputs. We observe that t-evolutions are slightly steeper in NLO calculations 
then those of LO. We can also calculate x-evolution of non-singlet structure functions at 
low-x from equation (10). But it involves complicated triple integrations and we keep it as 
our subsequent work. 

In figure 2 we plot T{tf and kT(t) against in the range 2 < < 59 GeV^ as 
required by our data used[15]. Though the explicit value of k is not necessary in calculating 
t-evolution of -^2^'^, yet we observe that for k= 0.027, errors become minimum and varies 
from 17.68% to -13.68% in the range 2 < < 50 GeV^. 
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Figure 1: t-evolutions of non-singlet structure functions F2^'^(a;, t) (solid lines) for the repre- 
sentative values of x given in the figures. Data points at lowest-Q^ values in the figures are 
taken as input to test NLO t-evoluation of non-singlet structure functions ^2^"^ from equa- 
tion (8). In the same figures we also plot the results of t-evolutions of non-singlet structure 
functions ^2^'^ (dashed lines) for LO from equation (9). 
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Figure 2: T{tf and kT(t) against in the range 2 < < 50 GeV^ for k=0.027. 
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